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Abstract 

We study convergence of the following discrete-time non-linear dynamical system: n agents 
are located in M. d and at every time step, each moves synchronously to the average location 
of all agents within a unit distance of it. This popularly studied system was introduced by 
Krause to model the dynamics of opinion formation and is often referred to as the Hegselmann- 
Krause model. We prove the first polynomial time bound for the convergence of this system 
in arbitrary dimensions. This improves on the bound of resulting from a more general 

theorem of Chazelle 0]. Also, we show a quadratic lower bound and improve the upper bound 
for one-dimensional systems to 0(n 3 ). 

1 Introduction 

Lately, there has been a surge of attention given to network-based dynamical systems, in which 
agents interact according to local rules via a dynamic communication graph [7]. Much of the 
previous work has focused on the exogenous case, where the communication topology is decoupled 
from the evolution of the system. We refer interested readers to [3j for a good overview of research 
in this area. In the more common, endogenous version, the communication graph changes over 
time according to the current states. The feedback loop between dynamics and topology creates 
considerable difficulties, and efforts have been underway to build an algorithmic calculus within the 
broad framework of influence systems |5]. This work investigates the complexity of Hegselmann- 
Krause systems (abbreviated as HK system from now on), a popular model of opinion dynamics 
that has proven highly influential over the years [U El El [TTJ [13] and stands as the archetype of a 
diffusive influence system [6j. In the d-dimensional version of the model, each agent has an opinion 
represented as a point in M. d . Two agents are neighbors if they are within unit distance from each 
other. At every time step, each agent moves synchronously to the mass center of its neighbors. 

HK systems are known to converge |XQ^ [T2l I16j . meaning that they eventually come to a full 
stop. The convergence time has been bounded by time and conjectured to be polynomial 
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It was shown to be 0(n 5 ) in the case d = 1 [13] • The contribution of this work is threefold: first, we 
prove that the convergence time is indeed polynomial in n, regardless of the dimension; second, we 
lower the one-dimensional bound from 0(n 5 ) to 0(n 3 ); third, we establish a quadratic lower bound, 
which improves the known bound of O(n) [Hj. We also consider noisy variants of the model. 

The bidirectionality of the system plays a crucial role in the proof, as it should. Indeed, it 
is known from [3] that allowing different radii and averaging weights for the agents can prevent 
the convergence of the communication graph. Our proofs are an elementary mix of geometric and 
algebraic techniques. Much of the current technology for HK systems centers around products of 
stochastic matrices. This work injects a geometric perspective that, we believe, will be necessary 
for further progress on the more difficult directional case. 



2 Preliminaries 

We formally define the discrete-time HK system in dimension d as follows. There are n agents. For 
every t G and for every i G [n], the position of agent i at time t is Xj(t) G M. d . The positions at 
t = are given and, thereafter, are updated synchronously according to the following rule: 

M t + l > = ^ j— (i) 

Here, || • || denotes the Euclidean norm. We say that agents i and j are neighbors at time t 
if ||xj(i) — Xj(i)|| ^ 1, and we denote the neighborhood at time t of agent i by Afi(t) = {j 
i and j neighbors at time t}. Note that if i G A/}(i), then j G ftfi(t). Also for a given x G 
we define the weight of x at time t to be Wt(x) = \{i : Xj(i) = x}|. Note that at any given t, 

^ ti^(x) ^ n for all x G ~K d and that the sum of all weights equals n. Also, let the weight of agent 

1 at time t denote tOt(xj(i)). The weight of an agent is monotonically non-decreasing with time. 
The system is said to have converged at time t if Xj(i + 1) = Xj(t) for all i G [n]. 

In the case of one dimension, one can observe that order is preserved: 

Proposition 2.1 (Lemma 2 in [H]). T/ie i^JC system in dimension 1 preserves order of positions. 
That is, for any i,j G [n], if Xi(0) ^ Xj(0), then Xi(t) ^ Xj(t) for all t > 0. □ 

For this reason, in one dimension, we can number the agents from 1 to n such that X\(t) ^ 
X2(t) ^ • • • ^ x n {t) for all t ^ 0. The one-dimensional system also has the following decomposability 
property: 

Proposition 2.2 (Proposition 2 in [3]). Suppose d = 1. // a£ time £, and /or some i £ [n], 
|xj+i(t) — Xi(t)\ > 1, then for all if > t also, \xi+\{t!)— Xi(t')\ > 1. So, tae system can be decomposed 
into two subsystems, one consisting of agents {1, . . . , i} and the other of agents {i + 1, . . . ,n}, each 
evolving independently after time t. □ 

Therefore, if at time t, an agent has no neighbor strictly to its left and no neighbor strictly to 
its right, it never moves subsequently, and we say that the agent has frozen at time t. 

This decomposability property is not true in higher dimensions. For instance, consider the 
two-dimensional HK system with agents at positions (0,0), (1,0) and (1/2, 1); the third agent does 
not neighbor any other agent at t = 0, but this is not so at t = 1. 



2 



3 Convergence in One Dimension 



Theorem 3.1. The HK system in 1 dimension converges within 0(n 3 ) time steps. 

Proof. Recall that we number the agents such that Xi(t) ^ Xi + ±{t) for all t and all i € [l,n — 1]. 
Suppose the system has not already converged at time t, and let £(t) denote the leftmost non-frozen 
agent, i.e., the least i ^ 1 such that {xi(t) : i G N(,{t)} / |a^(i)}. Note that at time t, agent £(t) 
must have at least one neighbor strictly to its right and no neighbor strictly to its left (because any 
neighbor strictly to the left would violate minimality of £). 
The following lemma is the heart of our proof. 

Lemma 3.2. For every t ^ 0, by time t + 2, agent £{t) increases in weight, or gets frozen, or moves 
to the right by at least 

Proof. Fix t, let i = £{t), and let r = min{j : Xj(t) > xg(t)} be the leftmost agent that is strictly 
right of agent £. Agent r is a neighbor of £ because £ has at least one neighbor strictly to its right. 
If Mi(t) = J\f r {t), agent £ at time t + 1 moves to the same location as agent r does and increases 
its weight. So, suppose otherwise. Since £ has no neighbor strictly to its left, r also does not have 
any neighbor strictly left of I. Hence, there must exist s G M r {t) \ftf((t) strictly to the right of r. 

We now show x r (t + 1) ^ xt(t) + \. Observe that x s (t) — xe(t) > 1. So, if S = x r (t) — xg(t) and 
k = \Af r (t)\, then agent r moves to the left by at most: 

5-{k-l)-{l-5)-l _ 5 1 
k k n 

If Xi(t + 1) ^ xt(t) + we are already done. Otherwise, x r (t + 1) — xg(t + 1) ^ By 
Proposition 2.2, agent £ still has no neighbor strictly to its left at time t+1. If x r {t+l)—xi{t+l) > 1, 
then agent £ gets frozen at time t + 1. Otherwise, ^ ^ x r (t + 1) — x^it + 1) ^ 1, and so: 

^ + 2) -^ + 1) ^2^TTl4 

□ 

We can assume x n (0) — x\(0) ^ n without loss of generality, because otherwise, by Proposi- 
tion 2.2, the system can be decomposed into independently evolving subsystems. So, x n {t) — x\(t) ^ 
n for all t. Now, apply Lemma 3.2 at t = 0, 2, 4, . . . as long as the system has not converged and 
£{t) exists. £{t) can increase in weight only at most n times. Also, because £(t) is non-decreasing 
with t, X£( t+ 2) {t + 2) ^ X£(f){t + 2) and, hence, the third case in Lemma 3.2 can occur only at most 
2n 3 times. Thus, after t > 2{n + 2n 3 ), Lemma 3.2 cannot be applied, and the system must have 
converged. □ 

3.1 Noisy neighborhoods 

In this subsection, we study an extension to a noisy version of the HK model. We consider the 
following system HK^ where r\ G (0, 1) is a parameter. There are n agents, and each agent has an 
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associated left-neighborhood parameter rji that is in the interval (0,r/). The positions at t = are 
given, and then, the positions are updated according to the following rule: 

a?i(t + l) = ^ j— (2) 

3':— l+»7i<asj(t)— Xi(t)<l 

One can interpret the r/j's as noise acting on the left side of each agent's neighborhood. We can 
prove the following extension of Theorem 3.1 to HK^. 

Theorem 3.3. For fixed r\ G (0, 1), the HK^ system converges within 0(n 3 ) time steps. 

Proof. The proof proceeds similarly to that of Theorem 3.1. We say that agent j is a neighbor of 
agent i at time t if — 1 + rji ^ Xj(t) — X{(t) ^ 1. Conceptually, the main complicating issue is that 
being neighbors is no longer a symmetric relation. Also, Proposition 2.1 is no longer true as order 
may not be preserved by the dynamics. 

Let £{t) denote the leftmost agent at time t that neighbors at least one agent positioned strictly 
to its right. Notice that this implies there is no k such that < x^m — Xk ^ 1, because agent 
k would neighbor £(t) and hence violate the minimality of £(£). We prove the following analog of 
Lemma 3.2. 

Lemma 3.4. For every t 0, by time t + 2, agent £{t) increases in weight, or gets frozen, or moves 
to the right by at least s-y. 

Proof. Let t = £(t) and let r be the leftmost neighbor of £. If x r (t) — x^{t) > 1 — rj r , then 
xi{t + 1) > 0-—Vr)/n ^ ( 1_? ?)/n, and we are already done. So, assume otherwise. 

Then, agent £ is a neighbor of agent r. Since £ has no neighbor strictly to its left within a 
distance of 1, r also has no neighbor strictly to the left of £. Now, we can use the same argument 
as in the proof of Lemma 3.2 to argue that either agent £ freezes at time t + 1 or at time t + 2, it 
moves to the right by at least ^2 • D 

The rest of the proof can be finished in exactly the same way as the proof of Theorem 3.1. □ 

By symmetry, Theorem 3.3 also applies when the right side of each agent's neighborhood is per- 
turbed and the left side is fixed. The case when both sides are perturbed remains open; convergence 
for such heterogeneous HK systems is conjectured [T5] , 



4 Convergence in Higher Dimensions 

In this section, we consider the HK system in d dimensions with n agents and d 2. We show that 
the convergence time is polynomial in both n and d. 

The proof follows from a sequence of lemmas. The first lemma asserts that there is some vector 
a such that taking the projection of agents on a does not bring two non-neighbors too close together. 



emma 4.1. For every t, there exists a unit vector a such that for any two agents i,j, we have 
Sl(n- 2 d- 1 ). 



_ Xj(t)-Xj(t) 



|[x»(*)-Xj(*)ll 
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Proof. Let dj = M x '^| Xi %n - We can view dj ,■ as a point on S d 1 , the unit ball in M. d . The 
set of points on S^ 1 with dot product with djj of absolute value at most 0(n^ 2 (i _1 ) has area 
n 2 dT((d-2)/2) • area °^ ^ 1 * s r(rf/2) • There are pairs i, j, and so, by a volume argument, 
there exists a point a G S^ -1 such that for all i and j, |dj j • a| = VL(n~ 2 d~ l ) . □ 

The next lemma analyzes the one-dimensional system formed by projecting onto a and shows 
a lower bound on the total movement in each step. 

Lemma 4.2. Assume that at time t, the system has not converged. One of the following cases 
happens. 

- Two agents move to the same location at time t + 1. 

- Some agent moves by at least Q(n~ d -1 ) 

Proof. Let a be a unit vector satisfying Lemma 4.1 and j be the agent with minimum Xj(t) ■ a 
among agents with neighbor at time t. Consider two cases. 

Case 1. j has a neighbor i such that (xj(i) — Xj(t)) • a = r2(n~ 3 d _1 ). Because (xfc(i) — Xj-(t)) • a ^ 
V/c, j must move by Q(n d ). 

Case 2. All neighbors i of j satisfy (xj(t) — Xj(t)) • a = 0(n _3 <i _1 ). If none of them has any 
neighbor that is not j or j's neighbors, then they all move to the same location at time t + 1. 
Otherwise, some neighbor i of j has a neighbor k that is not a neighbor of j. Then (xfc(i) — Xj-(t))-a = 
nin^d- 1 ). We have 

Mt + 1) • a > x*(t) • a + ^ Xi(t) + n^-aj-i) 



Thus i moves by at least Q(n d ). □ 

We now use the following special case of a result from [T7] to show the existence of a potential 
function which is strictly decreasing as long as some agent moves. The proof is included for 
completeness. 



/||2 

is 



Theorem 4.3 ([T7]. Theorem 2). Lei /| • : R d x M d -> R be defined as follows: //j(z,z') = ||z — z' 
j/ ||xi(t) -x 3 -(t)|| < 1 and flj(z,z') = 1 otherwise. Let V(t) = Ei Je[n] /^(x^'x^i)). T/ien F 
non-increasing along the trajectory of the system. In fact: 

V(t)-V(t + 1)> Yl /^•(xi(t + l)-x i (t),x j (t)-x i (t + l)) 
^^4||x,(t + l)-x^)|| 2 

i 

Proo/. First, notice that /f+^x^t + l),x,-(t + 1)) < /^(x^t + l),Xj(t + 1)) Vi, j. Next, we will 
show 
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M i,ie[n]:/*^l 

Let F : (R d ) n x (R d ) n ^ M be defined as F(x,y) = Eij^^i /*j( x i> Yi)- Then x + 1) = 
argmin y i ? (x(t), y). We can also view F as 2 matrices j4, i3 and -F(x, y) = x T Ax + y T Ay — 2~x T By 
where A,B are symmetric and A is positive semidefinite. We have x(i + 1) = ^4 _1 i?x(t). 

F(x,x) -i^yT^x^^Sx) = 2x T (BA~ 1 (A — B)A~ l B + A — B)x 

= 2x T (BA- 1 B - BA^BA^B + A- B)x 
= 2x T (/ - BA~ l )(A + B)(I- A _1 5)x 
= F((J - A _1 B)x, -(/ - A _1 5)x) 



□ 



Theorem 4.4. T/ie system converges in poly(n, d) time. 



Proof. There are at most n time steps where two agents move to the same place. In all other 
time steps, by Lemma 4.2, some agent moves by at least f2(n -4 <i -1 ). We have V(0) ^ n 2 and 
V(t) ^ Vi. Therefore, by Theorem 4.3, the number of time steps before the system converges is 
0(n w d 2 ). □ 



5 Lower bound 

In this section, we give an instance of an HK system that requires VL(n 2 ) steps to converge. Our 
example is in two dimensions, and we know of no example that takes longer time even when the 
number of dimensions is large. The previous best lower bound on the convergence time [14] was 
Q(n), which is achieved by n points on a line with unit distance spacing between each pair of 
consecutive points. This is still the worst example we know of in one dimension. 

Our lower bound is achieved by the following system. There are n agents in the system and 
initially, they are located at vertices of a regular n-gon with side length l\ = 1. Label the agents 
clockwise around the polygon from to n — 1. Let O be the center of the polygon and Ai be the 
location of agent i. For notational convenience, all index computation is done modulo n (so if i = 
then Ai_i is A n ^\). 

Theorem 5.1. The system described above requires at least fi(n 2 ) steps to converge. 

Proof. We will prove by induction that for all t ^ n 2 /28, A^s are vertices of a regular n-gon with 
side length l t > (1 - 

By the initial state of the system, the invariant holds for t = 1. Assume that it holds before 
step t = k, we will show it holds before step t = k + 1 ^ n 2 /28. 

We analyze the movement of agent i from Ai to A i in one step (see Fig. 1). Note that 
ZA{OAi + \ = 2 ^-. For any j £ {i — 1, i,i + 1}, we have 
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Figure 1: In one time step, agent i moves from Ai to A' { , the centroid of A{, A^\, Ai + i 
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Thus, for j {i — + 1}, agents i and j are not neighbors. Ai moves to the centroid of 
Ai, Ai_i, Ai + \ and by symmetry, the locations of the agents are vertices of a new regular n-gon 
centered at O and with a smaller side length h+i- We have 
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Thus, the system requires at least n 2 /28 steps to converge. □ 
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6 Discussion 



In this paper, we analyzed the convergence rate of the homogeneous HK system in arbitrary di- 
mensions. The system is shown to converge in polynomial time, but can take at least a quadratic 
number of steps in the worst case. Getting tight bounds on the convergence time of the system, 
even in just one dimension, remains an interesting open problem. 

A particularly interesting challenge is to analyze the heterogeneous version of the HK system, 
i.e. when the neighborhood radii of all the agents are not necessarily the same |15 |. New ideas 
are needed to understand the behavior of this system in particular, and directional systems in 
general. Our analysis of a noisy variant of homogeneous HK system is a step towards studying 
more complicated directional systems. 

Beyond convergence rate, the behavior of the homogeneous HK system is still full of mysteries. 
Most notable perhaps is the 2R conjecture [2]: when agents are drawn uniformly at random on 
an interval, they converge to clusters at distance close to twice the minimum possible inter-cluster 
distance. Resolving this conjecture remains well beyond our understanding of the system. 
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